Quantum phase slips in the presence of finite-range disorder 
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To study the effect of disorder on quantum phase slips (QPS) in superconducting wires, we 
consider the plasmon-only model where disorder can be incorporated into a first-principles instanton 
calculation. We consider weak but general finite-range disorder and compute the formfactor in the 
QPS rate associated with momentum transfer. We find that the system maps onto dissipative 
quantum mechanics, with the dissipative coefficient controlled by the wave (plasmon) impedance Z 
of the wire and with a superconductor-insulator transition at Z = 6.5 kOhm. We speculate that 
the system will remain in this universality class after resistive effects at the QPS core are taken into 
account. 

PACS numbers: 74.78.Na, 73.21.Hb 
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The possibility that quantum fluctuations destroy su- 
perconductivity in thin wires has attracted attention of 
both experimentalists and theorists for a long time.JSim- 
ilarly to Little's analysis of thermal fluctuations, El one 
concludes that the requisite quantum fluctuation should 
be sufficiently large, so as to allow the Ginzburg-Landau 
(GL) order parameter t/jql to vanish at the core, and the 
phase of ipGL to unwind. Such fluctuations are known as 
quantum phase slips (QPS) 13 On the experimental side, 
there has been a bit of controversy over precisely how su- 
perconductivity disappears in thin wires at low temper- 
atures. Some experimentsj-see a sharp superconducting- 
insulator transition (SIT), til while others do not.Q 

Superconductivity takes place when it is difficult to 
transfer momentum from the moving condensate. In 
one dimension (ID), a QPS unwinds a large momentum 
P ~ nn s , where n s is the linear superconducting elec- 
tron density, and this momentum has to go somewhere. 
In uniform superfluids (such as a cold Bose gas in a ring 
trap), the requirement of momentum transfer constitutes 
a major bottleneck for QPS.B In superconducting wires, 
on the other hand, there are some easily identifiable sinks 
of momentum. The most obvious, and as far as we know 
the only one that has been considered in the literature, 
is the normal electrons at the QPS core, which in turn 
transfer momentum to the disorder potential. It appears 
that in the existing theory of QPS in wirestl this process 
is assumed to be 100% effective, so that no trace of mo- 
mentum conservation is left in the QPS rate. One should 
keep in mind, though, that this result is obtained using 
instantons of a disorder-averaged theory. 

In this work, we analyze effects of suppressed momen- 
tum transfer explicitly, within a simple model where the 
QPS rate can be found from a first-principles instanton 
calculation — by first obtaining the rate for a given disor- 
der configuration and then averaging over disorder. We 
consider the general case of weak but finite-range disor- 
der. Our essential simplifying assumption is that elec- 
trons in the core are not effective in transferring momen- 
tum to the lattice, so thai the transfer takes place via 
the gapless plasmon mode.13 In this limit (applicability of 



which is further discussed below), the rate can be com- 
puted within a plasmon-only effective theory [Eq. (|l|)]. 

Under these assumptions we find that the transition 
point is determined by the wave (plasmon) impedance 
of the wire Z. We find that the system is in the .-uni- 
versality class of the dissipative quantum mechanics^ (as 
opposed to the XY universality class found in Ref. ||) and 
identify a SIT at Z = n/2e 2 = 6.5 kO. The Ohmic resis- 
tivity of the wire at the SIT, for weak disorder (the only 
case considered here), is small (much smaller than the 
normal-state resistivity) and has a non-universal value 
that depends on both the strength and correlation length 
of disorder. 

These results apply in the limit when the normal re- 
sistance i? C orc of the QPS core is effectively infinite. The 
effect of a finite i? CO rc can be understood as follows. Plas- 
mons produced by a QPS can be viewed as charge fluctu- 
ations in equivalent transmission lines, one such line on 
each side of the QPS core. A finite i? C ore will shunt the 
charge separation at the QPS core, thereby reducing the 
plasmon emission. This picture of two transmission lines 
shunted by i? cor e suggests that the universality class will 
remain the same even when dissipation is caused mostly 
by a finite i? CO rc (which may very well be the case for 
existing experimental samples). The SIT will now be 
controlled by the total impedance formed by Z and i? C orc 
connected in parallel and thus occur across a straight line 
in the (l/Z, l/R COTe ) plane. 

Our results rely on a certain amount of impedance 
matching at the wire's ends. We assume that plasmons 
can leave the wire and go into the leads. This inhibits 
quantization of plasmon modes and translates, techni- 
cally, into the possibility to consider the temperature T 
and the wire length L w as independent infrared parame- 
ters. In the limit of short-range disorder, our results can 
be compared to those obtained by Luttinger-liquid meth- 
ods in Ref. || and reproduced in the instanton approach 
in Ref. [l0| Because we use the scaling T — > with L w 
fixed (as opposed to T a L^ 1 used by those authors), 
we obtain a different value of the critical coupling. On 
the other hand, our value of the critical coupling coin- 
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cides with that obtained in the model where an effective 
resistor is connected to the wire's ends,M provided we 
substitute Z for the resistance. 

We start with the purely bosonic Euclidean lagrangian 
density 
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L E = ^d T ^ + ^\d x ^\ 2 + ^\ 4 



(»+vm\ (i) 



essentially the ID Gross-Pitaevskii model in the presence 
of the disorder potential V = V{x). The field tp is the 
"field of Cooper pairs" describing fluctuations of super- 
conducting density and phase. Thus, A^t^m^/M = n s 
(to* is the effective electron mass) is the density of su- 
perconducting electrons, while the coupling constant in 
Eq. (0) is g — 4e 2 /C, where C is the wire capacitance per 
unit length. The effective theory (0) holds only at spatial 
scales larger the size of a Cooper pair, i.e., the GL coher- 
ence length £. Thus, the potential V(x) is coarse-grained 
at the scale £ in the x-direction and at the scale of the 
wire thickness in the transverse directions. 

We assume the disorder to be Gaussian with the corre- 
lators (V) = 0, (V(x)V(x')) = V$f{x-x'). The correla- 
tion function f(x) is normalized so that /(0) = 1, and Vq 
is the r.m.s. disorder amplitude. The disorder correlation 
length is meaning that for x > I, f(x) — > 0. 

A major role in determining the QPS rate is played by 
interactions between QPS at different locations. These 
interactions are determined by regions outside the QPS 
cores, and can be accounted for in a phase-only model. 
We define ip = (p + SpY^e 1 ® and, assuming weak disor- 
der, expand Eq. (|l|) in powers of the small density fluctu- 
ation 5p around the stationary point with a given phase 
gradient 9' — 9'. This phase gradient takes into account 
the biasing current / = 2e9'p/M. At the classical level 
the stationary point is characterized by the local mini- 
mum p — p± of the effective potential 



u (p) = T,P 2 - PP - 



Ml 2 
8e 2 p' 



(2) 



The minimum exists below the critical current, I < I c , 
where 3[MI 2 g 2 /(ie 2 )} 1 / 3 = 2p. 

Integrating out the density fluctuations Sp, we obtain 
the Euclidean lagrangian density for the phasa-fluctua- 
tions 9\ =9 — 9, describing gapless plasmonsQ propa- 
gating with speed c = {p+g/M) 1 / 2 on the background 
moving with superfluid velocity u = I /(2e p-,,): 

L E = ipA + X -B r 9 x K~ x D r 9x + ^{9[+ 9') 2 , (3) 

where D T 9\ = 9\ —iu9' x —iV is the covariant time deriva- 
tive in the moving reference frame and K = g (1 — r 2 V 2 ) 
is the differential operator with the screening length 
r s = (4M p+g)^ 1 / 2 . For realistic values of the parame- 
ters, this screening length is much smaller than the GL 
coherence length £. Our starting point ([!]) is already 
coarse-grained at scale £; in what follows we set K = g. 



We can now find the exponential factor in the QPS rate 
by computing the action of a suitable classical configura- 
tion. We begin with the case of strictly zero temperature, 
T = 0. The leading effect is due to a single phase slip- 
antislip pair, or equivalently a vortex-antivortex pair in 
the (x, t) plane. Away from the vortex cores 



x - to_(t - t ) 



w_(t - r' a ) 



. (4) 

where v^i = Co T u are the up(down)-stream velocities. 
For u — (i.e., v+ — f-), this is the configuration famil- 
iar from the studies of the planar XY model£3 

Integrating the Euclidean lagrangian density (||) with 
the configuration (^) over x and t, we obtain the corre- 
sponding classical action 

Se = Si + S2 + Sdis, (5) 
where the combination of uniform linear in 9\ terms, 



Si = iPAx - EAtq, 



(6) 



with Axo = x — x' Q , At = t — Tq, accounts for the 
Berry phase of each QPS [P = lirp* is the momentum 
released by unwinding the supercurrent] and the released 
energy E = 2-xp+u = 2nl/(2e), 

irc r (Ai + iv + AT )(Ax - to_At )-, 
S 2 = ln[ 5 7 

9 r c 

is the plasmon-mediated interaction between the phase 
slips coming from the terms quadratic in 9\ , and 
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a{i-u 2 /c 2 ) 



dxV{x) 



(8) 



is the effect of the disorder. In Eq. (0), we have used the 
QPS core size rc (which is not determined by the present 
theory) as the short-distance cutoff. 

These expressions illustrate the effect of the super- 
fluid velocity u on the QPS action and can be useful, 
for instance, in weakly non-ideal Bose gases, where u 
may in principle approach cq- Thus, for example, factor 
(1 — u 2 /cq) _1 = (v^+vZ 1 ) c /2 in Eq. (||) arises because 
the up(down)-stream plasmons spend longer (shorter) 
time at the place with the given density. In thin su- 
perconductors, however, we typically have u <C vp cq. 
In the following we neglect u/cq <C 1 , to get 

S E = Sl + ---ln^--^-i P'dxVix), (9) 
9 r c g J x ' o 

where Ar = (Axq + CqAt 2 ) 1 / 2 . We observe that the last 
term here can be interpreted as the modification of Si due 
to the local correction to the density 5p+(x) = —V(x)/g 
caused by disorder. We will assume that V(x) incorpo- 
rates all mechanisms leading to linear density inhomo- 
geneity: non-uniform wire cross-section, magnetic impu- 
rities, etc. The weak-disorder approximation is applica- 
ble for 1 5 p^ I <C p±, which gives the dimensionless measure 
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of the disorder strength, 

a = Vo/Mcl = Vo/gp± 



analytic continuation. We obtain 



(10) 



For a given disorder configuration, the QPS rate (per 
unit length) can be found as the imaginary part of the 
partition sum of the pair, 

U=^4r, ^Sm/^ [dAT e- s °, (11) 
9 r c J -^w J 

where the dimensionful prefactor with the coefficient 
C\ ~ 1 incorporates the fluctuation determinant and the 
Jacobian of transformation to the collective coordinates. 
The integration over xq, x' is extended over the length 
L w of the wire. For large enough L w the rate ( |TT| ) is 
self-averaging with respect to disorder. In this case, the 
effective action is obtained by disorder averaging, 



27tc Ar 1 2 2 

S c s = Si H In h -P a A(x 

g r c 2 

where the function 



A(<c) = 2 dx'(x- x') f(x') 
Jo 



4), (12) 



(13) 



is proportional to the average square of the disorder- 
induced phase between the slip and antislip. [The dis- 
order correlator f(x) is defined after Eq. (jl|).] For finite- 
range disorder, A(x) has a diffusive form at large dis- 
tances, \x\ >• I, 



A(x) = 2h\x\, l x = / dxf(x). 



(14) 



For single-scale disorder, lj ~ I. 

According to Eqs. (12), dl4|), disorder binds together 
the x coordinates of the phase slip and anti-slip within a 
pair. For the case of short-range disorder, this effect was 
noted previously in Ref. ^ and was taken into account by 
using a sharp <5-function, 5(xq — x' ), in the partition sum. 
We see, however, that in general there is a finite charac- 
teristic separation Axo within the pair. This separation 
depends both on the range and on the strength of the 
disorder potential. As the strength of disorder increases, 
and the pairs become more tightly bound in the x di- 
rection, the wire effectively becomes more and more like 
a dissipative quantum-mechanical system with a loga- 
rithmic interaction, In \tq — Tq|, between the instantons — 
a spatially extended version, of the familiar resistively- 
shunted Josephson junction.El 

We also note that the coefficient in front of the loga- 
rithm in Eq. (|l^) is proportional to the inverse impedance 
of the wire viewed as two transmission lines attached to 
the phase-slip region. This is consistent with the map- 
ping onto dissipative quantum mechanics, since plasmons 
are the only source of dissipation in our present model. 

Turning to the time integration in Eq. ([IT]) , we observe 
that while the integral is formally divergent, the corre- 
sponding imaginary part is finite and can be found by an 
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r(z/ + i/2) ' 

(15) 

where J„(z) is the Bessel function and the dimensionless 
index v is inversely proportional to the interaction con- 
stant, v + 1/2 = nco/g (= K^ 1 in notations of Ref. ||). 
This explicit expression for the QPS rate (at T = 0) 
in the plasmon-only theory in the presence of weak but 
finite-ranged disorder is the main result of this work. 
Although Eq. ([l5]) was obtained assuming that disor- 
der induces only weak fluctuations of the density (as 
expressed by the weak-disorder condition a -C 1), it is 
non-perturbative with respect to the phase fluctuations 
and accounts for multiple scattering to all orders. 

In the absence of disorder [a = 0] the integral in 
Eq. (|l5|) is zero for any E < cqP, consistent with the 
Galilean invariance of the T = state.B The opposite 
case is when the convergence of the integral (|l^) is dom- 
inated by disorder. In this case, the Bessel function can 
be replaced by the first term of its small-argument ex- 
pansion, and the QPS rate becomes 



g 2 r A c V 2c / J h 



iPx- 



2 A(x) 



(16) 

with C 2 = 7T 3 / 2 Ci/[r(z/ + \)T{v + 1/2)]. This limit cor- 
responds to the entire momentum P being absorbed by 
the disorder, with no momentum carried away by plas- 
mons. For weak disorder, where Eq. ( jl6| ) is applicable, 
the dimensionless formfactor is small. 

At nonzero temperature, the simplest case is T ^ 
gE/ircg. Then, the energy released by unwinding the 
supercurrent is insignificant, and instead of a single pair 
(^) we can use a periodic chain of such pairs with pe- 
riod 8 = l/T. Properties of this chain are described in 
Ref. Q An especially simple result applies in the disorder- 
dominated regime, when the spatial separation in each 
pair is small, Aa;o *C cq/ttT: for an estimate, it is suffi- 
cient to make the replacement E — » ircoT/g in Eq. (p^). 
Then, at v ~ 1, the rate of thermally-assisted QPS per 
unit length is 



r 6 G \ g J 

The voltage across the wire is 

2?r ttI 
V = — TZ T L W smh — — . 
e 2eJ 

We see that resistance becomes T-independent at v = 
1/2, which defines the SIT point. The dimensionless mea- 
sure of QPS pair density near v = 1/2 is 



(17) 



(18) 



IZtAxqAtq 



hAx /TrcY"- 1 



•c 



(t 1 )' 



(19) 



where Atq — 1/2T is the characteristic size of a pair 
in the r direction. Thus, the density is small in the in- 
frared at any v > 1/2, and for sufficiently small A A even 
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at v = 1/2. At v < 1/2, pairs proliferate, resulting in 
the insulating behavior. Writing the impedance of the 
wire as Z = 2 (L/C) 1 / 2 , where L = m*/e 2 n s is the "ki- 
netic" inductance per unit length, and the factor of 2 
corresponds to two transmission lines (one on each side 
of the phase slip), we see that v = 1/2 is equivalent to 
Z = 7r/2e 2 = 6.5 kfi. 

We now turn to discussion of the formfactor The 
important parameter here is a 2 Pl. We begin with the 
case a 2 PI <C 1. If the stronger condition aPl < 1 is 
also satisfied, we can (upon integrating by parts twice) 
expand the integrand in ( [l6| ) to the linear order in 
a 2 . This corresponds to the entire momentum P be- 
ing absorbed in a single scattering event. Then, = 
(a 2 /h) J e lPx f(x)dx. Thus, for PI <C 1 the formfactor is 
universal, A& = 2a 2 , while for I > P" 1 it depends on de- 
tails of the correlations. For PI 3> 1, it is determined by 
the ordinate Iq of the singularity of f(x) closest to the real 
axis and scales exponentially, A^ oc e~ pl °. If aPl > 1 
(but still a 2 PI <C 1), effects of multiple scattering result 
in a correction ~ a 2 P 2 l 2 in the exponent. 

In the opposite limit, a 2 Pl ;» 1, it is the disorder 
that determines convergence of the integral (|l6|). The 
integral converges at x ~ (a 2 P) _1 <C I, where we can 
replace A(x) by its short-distance form, A(x) = x 2 . 
Then, the integral becomes Gaussian and gives A& — 
{2-K) 1 ' 2 (aPl\)~ l e~ 1 /( 2a '. In this case the process of mo- 
mentum transfer to disorder is clearly a result of a large 
number of scattering events. 

The nature of the crossovers between different scat- 
tering regimes as the disorder correlation length I is in- 
creased can be understood in a model with correlation 
function f(x) = (1 + x 2 /l 2 )- 3 / 2 . Integration in Eq. |[|) 
gives A(x) = 2/ 2 [(l + x 2 /l 2 ) 1/2 - 1], and the coordinate 
integration in Eq. (|l5|) can be performed explicitly, 



A d = 



2a 2 Pl 



[1 + (a 2 Pi) 2 ] 1/2 



^111 + (a 2 PI) 2 ] 1 ' 2 ). 



For PI <C 1, we can use K\(z) = z 1 + O(zlnz) which 
gives the short-range limit (with some logarithmic cor- 



rections due to the power-law tail of the correlation func- 
tion), while for PI ^> 1, the large- argument asymptote 
K\(z) w (n/2z) 1 / 2 e~ z can be used to restore the other 
discussed scattering regimes. 

Our discussion of A d has so far assumed that the 
scattering is characterized by a single distance scale, I, 
which may not necessarily be the case. For example, a 
model correlation function f(x) — e~\ x \l l has a singular- 
ity (derivative discontinuity) at the real axis, which gives 
distance Iq = 0, and the exact integration in Eq. ( |l6| ) 
shows that A d is independent of P for the entire range 
aPl <C 1. Conversely, for a correlation function with the 
power-law long-distance tail, f(x) oc x~ m , < m < 1, 
the distance l\ is infinite, and instead of Eq. ( |i~4| ) one 
obtains superdiffusive form 



A(x) 



(-)' 



x ^> h 



2 r TO 

where is a convenient distance scale. Then, pertur- 
bation theory breaks down already for arbitrarily weak 
disorder; we get ln(A d ) oc ~[a 2 (Pl lf ) m }- m/{1 ~ m) ■ 

In conclusion, the plasmon-only theory provides a use- 
ful laboratory for studying the effect of disorder on QPS 
from first principles. We have computed the QPS rate for 
general finite-range weak disorder; as seen from Eq. (|l|), 
disorder drives the system into the universality class of 
dissipative quantum mechanics, with the dissipative coef- 
ficient determined by the wave impedance Z of the wire. 
From an extension of this result to finite temperatures, 
Eq. ( |l7| ) , we have found a superconductor-insulator tran- 
sition at Z = ir/2e 2 . We see no reason why the universal- 
ity class should change when resistive effects at the QPS 
core are taken into account. Rather, we expect, that it 
will remain the same, but the dissipative coefficient will 
now be determined by the total impedance, including 
both plasmon and resistive effects. 
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